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Abstract 

We consider the problem of obtaining integral representation of feedback operators for damped hyperbolic 
control systems. We show that for the wave equation with Kelvin- Voigt damping and non-compact input 
operator, the feedback gain operator is Hilbert-Schmidt. This result is then used to provide an explicit 
integral representation for the feedback operator in terms of functional gains. Numerical results are given 
to illustrate the role that damping plays in the smoothness of these gains. 
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1 Introduction and Motivation 

The purpose of this article is to extend the representation theorem in [1] and [7] to certain classes of damped 
hyperbolic systems. The original motivation for our study of hyperbolic systems comes from the work by 
Lupi, Chun, and Turner [8]. The approach in [8] is interesting because they make no prior assumptions 
regarding the form of the controls and actuators so that the gains operators produced by an optimal design 
could be used to make decisions about where actuators and sensors are best placed. In particular, in [8] it 
was assumed that the input operator was the identity and the elastic system was not damped. By solving the 
LQR problem with the input operator equal to the identity, one can gain insight into the type and location 
of practical distributed controllers for structural control. This insight comes from explicit knowledge of the 
kernels (so called functional gains) that describe the integral representations of feedback gain operators. 
Even with no damping the LQR problem has a solution since the input operator is the identity (the system 
is exactly controllable). However, as we see below the problem with no damping is extremely complex. Basic 
questions concerning the existence and smoothness of functional gains remain open and yet these issues are 
important in the applications proposed in [8]. Consequently, as a first step we take the middle ground and 
consider damped systems with distributed control. 

Except for the obvious cases with bounded input operator and a finite number of controllers, the problem 
of obtaining explicit representations of feedback laws is more complex than one might first imagine. In the 
most general case, this problem is equivalent to the problem that led Grothendieck to develop the theory of 
topological tensor products and nuclear spaces. This theory led to the famous Schwartz Kernel Theorem. 
However, the fact that the operators of interest often arise as solutions to Riccati equations can be exploited 
to yield reasonable results. The representation problem for a parabolic problem with unbounded control 
operator was first considered in [1] and [7]. It was shown that as long as the input operator is bounded 
relative to the open loop dynamic operator, the solution to the algebraic Riccati equation is Hilbert-Schmidt. 
This fact was then exploited to show that the resulting feedback operator had an integral representation. 
For the hyperbolic case considered here, these types of results are more intricate and highly dependent on 
the type of damping. 

In the Section 2 we present a numerical example involving the control of a hybrid cable-mass system. 
This example is used to motivate the model problem and to demonstrate that distributed controllers can 
enhance disturbance attenuation. We then concentrate on the ID wave equation with Kelvin-Voigt damping. 
Although similar results can be obtained for other damping models, we present the basic theorem for this 
model in order to keep this paper moderately short. However, we illustrate the ideas and difficulties for 
other damping models with numerical examples. Finally we close with a few comments about future work 
and other open problems. 

2 Control for a Cable-Mass System 

The following system was proposed by Nayfeh, Nayfeh and Mook [11] as a simple example of a nonlinear 
distributed parameter system with the property that many standard discretized lumped models failed to 
capture the essential nonlinear behavior of the dynamic system governed by the partial differential equation. 
This system was also considered in [2] where it was used as a test model for MINMAX control for distur- 
bance attenuation. However, in [2] there was only one controller and, although disturbance attenuation 
was achieved, we shall see below adding a distributed controller can improve performance. The MINMAX 
approach provides a “robust state feedback control law” which is less sensitive to disturbances and certain 
unmodeled dynamics than is the LQR design. The idea is to obtain a representation of the control law 
and then use approximation theory to compute finite dimensional suboptimal controllers. These suboptimal 
controllers were used to attack the problem of designing reduced order state estimators. We shall limit our 
discussion here to the full state feedback problem. 

Consider the hybrid nonlinear distributed parameter system described by a vibrating cable held fixed at 
one end and with a mass attached at the other end. The mass is suspended by a spring which has nonlinear 
stiffening terms and is forced by a disturbance (see Figure 2.1). The equations for the hybrid system are [2] 
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The initial conditions are given by 


w(t , 0) = 0. 


(2.3) 


w(0,s) = w 0 , —w(0,s) = w 1 . 


(2.4) 


Here, u>(f,s) represents the displacement of the cable at time t, and position s, w{t,l) represents the 
position of the mass at time f, p and m are the densities of the cable and mass respectively, r is the tension 
in the cable, and 7 is a damping coefficient. The alphas are coefficients describing the nonlinear effects of the 
spring. The term jj(<) is viewed as a disturbance and ui(f,s) and «2 (t) are control inputs. For the moment 
we assume Kelvin- Voigt (internal) damping in the cable. 

The problem is hybrid in that the system is described by a linear partial differential equation (the wave 
equation) coupled through the boundary condition to a low order nonlinear ordinary differential equation 
(Duffing’s equation). In [2] it was assumed that the control acted exclusively on the mass (i.e. that ui(f,x) 
= 0 ). 

This model is often first written as a second order system in a Hilbert space H of the form 


y(t) + D 0 y(t) + Aoy(t) + F 0 (y(t)) = Bou(t). 

For the cable-mass problem considered here, H = 1) x H 1 and y(t) = [u>(<, •)> w(t, i)] T - 



I 

w(t,s) 


(2.5) 


Figure 2.1: Cable-Mass System 

This formal system has the advantage that it has the same appearance as the finite dimensional case 
and in order to address viscous and “structural” damping one merely replaces Dq = A 0 with Do = I and 
Dq = [Ao] l f 2 , respectively. We note however, that it is more consistent with physics to write the system in 
the form 

y{t) + S* (Sy(t) + iTy{t)) + F 0 (y(t)) = B 0 u(t) (2.6) 

where S — T — [y4.o]^^- Observe that Ao “ ^4q 0 and so S* = S and S* S — ST — Ao . Hence, (2.6) is 

formally obtained by factoring [. A 0 ] 1/2 out of the expression D 0 y(t) + A 0 y(t) in (2.5). Note also that (2.6) is 
of a form that allows for structural damping where S = [^o] 1/2 and T - 7, as well as for viscous damping 
where S — [Ao] 1 / 2 and T = S~ l = [Ao]~ 1 / 2 . In addition, by writing the system in the second order form 
(2.6), one captures a form that comes from balance laws and at the same time sets the stage for a simple 
formulation of the problem in first order state space form. 
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The system governed by equations (2.1) - (2.4) can be written as a dynamical system in an appropriate 
(infinite dimensional) state space. Although there are several equivalent formulations for this problem, we 
shall write the governing equations as the first order system 


x(t) = Ax(t) -j- F(x(t)) + Bu(t) + Dr)(t ), x(0) = x 0 (2.7) 

where at time t the state x(t) = [y(f ), y{t)] T lies in the Hilbert space X = Hi x IR x L 2 x 1R. Here, Hi is the 
subspace of the Sobolev space H l = H l [ 0,/] defined by = {w € i / 1 : tu(0) = 0}, and L 2 is the standard 
Lebesgue space of square integrable functions. The control u( 2 ) lies in the control space U = L 2 x IR. Here 
the inner product in X is 


{[^(•)> £>® 0 > p ] T » [*(•). = T [ w'(x)w'(x)dx + a & 

Jo 

+p I v f (x)v *{x)dx + mpp. 

Jo 


(2.8) 


It is important to precisely define the system operators and their domains in order to obtain correct rep- 
resentations of the feedback operators that will be used to control the system. Let Si denote the “evaluation 
operator” defined on H l [0,£] by Si(<f>(-)) = </>(£) and define the linear operator A on the domain V(A) C X 
by 


V(A)=^x = [vj ) ( t ,v,ii\ T € X :w,v e Hi, 1 “^™+ 

w(() =t,v(t) = /i}, 


and 


[ d f r d 7 (f 1 . f r (i 7 d 1 a\ . 

Ax= v,fi, — < w+ v > ——w + — — w > ^ 

as 1 p a$ p ds J [m ds mas J m 

The control input operator B and the disturbance operator D are defined by 


B — [0, Jff] T and Dr} — 
respectively. The nonlinear operator F is defined on X by 


n T 


0 , 0 , 0 ,— 7 } 
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0, 0, 0, K] 3 1 = [0, Fo(y)] 5 
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(2.9) 


(2.10) 


( 2 . 11 ) 


( 2 . 12 ) 


Observe that the input operator is the same as that used in [ 8 ] in their analysis. 

As noted above, this problem with Ui(t,s) = 0 was considered in [2] where MINMAX control was used 
to design a low order dynamics control law. This law was based on two outputs (position and velocity of 
the mass) and resulted in a practical low order design. For this note, we shall consider only the full state 
feedback problem. However, we allow for distributed control through ui(f,s). The simplest approach is to 
linearize the system, use MINMAX design to obtain a feedback operator and apply this law to the full plant. 

The linearized system has the form 


d 


dt 


( t ) 


Ax(t) + Bu(t) + Dr)(t), x(0) = xq. 


(2.13) 


For this problem, with Kelvin- Voigt damping, one can apply the MINMAX theory in [10] to obtain a feedback 
law of the form 


u(i) = 


«i (t,s) 

«2 (f) 


= -K,x(t) 


(2.14) 


where for 0 > 0 the gain operator K$ : X — > H is given by 


K e = 


(2.15) 
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(2.16) 


and Pg satisfies the algebraic Riccati equation 

A* P + PA - P[BB * -9 2 M]P + Q = 0. 

Here M = M* > 0, Q = <?* > 0 and the system (2.16) holds in the weak sense (see [10]). When 0 = 0 
one has the LQR design. However, when 6 > 0 the corresponding MINMAX controller provides additional 
disturbance attenuation (see [2]). 

Figure 2.2 clearly illustrates the difference between the performance achieved in [2] and what one can 
achieve with additional distributed control. For this example, we use the same parameters and finite element 
scheme found in [2]. 
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When there is control on the mass only, the maximum value of 6 yielding a feedback law of the form (2.15) 
- (2.16) was 9 = 1.7; this behavior is shown by the dashed line. If in addition, one allows distributed control 
(at 6 = 1.7) there is a loss in attenuation (but possibly better performance) as shown by the dotted line. 
However, by allowing control on the cable, the value of theta can be increased to 9 — 2.5, greatly increasing 
disturbance attenuation and improving performance as shown by the solid line in Figure 2.2. 



Figure 2.2: Displacement of the Mass under Three Types of Control. 


This example clearly indicates the benefit of distributed control and MINMAX design. Moreover, as 
shown in [2] if one has explicit representation of the feedback operator K e > then this can be used to design 
practical low order dynamic controllers (nonlinear compensators). In particular, it follows that since 




Kl 3 

K} a 

K? 

Kf 

K? 

Kf 


(2.17) 


is bounded from X to H , the control on the mass has the form 

r l Q 

u 2 (t) = — J rk s (z)—w(t,z)dz — aik a (t,£) 

- J pk v (z)-^w(t,z)dz - 


(2.18) 
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It is tempting to assume that one also has a representation of the form 

ft fS 

«i(*><) = ~ J Th s (C,z)—w(t,z)dz - cnh a (C)(t,e) 


ft ^ 

~ / pK{C,,z)—w{t,z)dz -mh m (C,) — w{t,t). 
Jo 


(2.19) 


dt 


dt 


In (2.18) k,(z) and k v (z) are functional gains corresponding to strain and velocity, respectively. Their 
existence and smoothness properties are assured by the Riesz Representation Theorem. For example, there 
is a k,(z) € Z-2(0,f) so that Kg 1 : H £ — + IR has the form 


f 1 J 

\ K e X ] 4>{-) = j Tk,(z) — <t>(z)dz. (2.20) 

However, on the surface all we know about K] 3 : 1,2(0, 1) — ► f-2(0, 1) is that it is bounded. As noted 
above, the desire to find representations of such operators as integrals led to Grothendieck’s work on nuclear 
spaces and the Schwartz Kernel Theorem. Recall that not all bounded linear operators on 1/2(0, 1) (even if 
self-adjoint) have integral representations, as illustrated by the identity operator. 

In [8] this issue was avoided by assuming a representation similar to (2.19) and then allowing generalized 
functions as kernels. This approach proved to be satisfactory for the one dimensional case considered therein, 
but does not apply to more general two and three dimensional hyperbolic problems. We present some results 
that lead to a representation of the form (2.19) when there is suitable damping in the system. Although this 
approach does not apply to the undamped case, the ideas can be extended to certain higher dimensional 
damped elastic systems. Moreover, we conjecture that the undamped problem in 3D systems will not yield 
a representation even of the type considered in [8]. 

We turn now to the simple ID wave equation in order to state precise results and to keep this article 
at a reasonable length. Moreover, we restrict our presentation to LQR design (9 = 0) since we can rely on 
existing literature to outline the results. Extensions to higher dimensional problems with 9 > 0 will appear 
in a future paper. 


3 The Wave Equation 


Consider the wave equation with Kelvin- Voigt damping. Damping plays a key role in the design of controllers 
for hyperbolic systems. To illustrate this role, consider the LQR control problem for the wave equation 
defined by 


I i* *(*■ s) = h x w{t ’ s) + s) + u(< ’ s)> 

0 < s < 1, 0 < t, 


(3.21) 


with boundary conditions 
and cost function 


w{t, 0) = 0, w{t, 1) = 0, 0 < t 


1 r) 

m = iL l 


di wit ’ s) 


+ M<.s)l 2 


dsdt. 


(3.22) 

(3.23) 


This problem is defined on the state space X = H q( 0, 1) x L 2 ( 0, 1). As above, care must be used to define 
the system operators. In particular, let 


and define A by 


V(A) = 


w 

V 


6 X : w, v € Hq(0, 1), w + jvE 


^(0,1)} 


A 


w 



»(•) 

V 


. ds 



(3.24) 


(3.25) 
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The control space is U = 1-2(0, 1) and the control operator is defined from U into X = Hq x L 2 by 


Here, 0 = 0 and Q - R = I X - The LQR problem has a solution (even if 7 = 0) given by 




w{t,-) 

&»(*.-) 


where K : Hl( 0, 1) x Z-2(0, 1) has the form 


K = B*P 


and P satisfies the weak form of the Riccati equation (ARE) given by 

( Px , Az)x + {Ax, Pz)x — (B m Px , B* Pz)u + ( Cx , Cz)y = 0, 
for all x,y in V{A). For this second order system, K takes the form 


or equivalently, 


where 


K = [ 0, 1l 2 ]P, 


K = [P 2 i, P 22 ] 


Pll P\2 

P21 P22 


(3.26) 


(3.27) 


(3.28) 


(3.29) 


(3.30) 


(3.31) 


(3.32) 


The operators P 2 1 : Hl(0, 1 ) -* L 2 (0, 1) and P 2 2 : L 2 { 0 , 1) — L 2 (0, 1) are bounded linear operators with 
p* - p 22 . The goal is to determine if there exist “nice” integral representations of these operators. 

The following result is well known and may be found in [3]. 

Lemma 3.1 The operator A generates an analytic semigroup S(t ) on X and there exist M > 0,w > 0 with 

||5(<)ll < Me-"*. 

Since A generates a stable analytic semigoup, we can apply Theorem 2.1 (page 36) in [9] to obtain the 
following regularity result for P. 

Lemma 3.2 There exists a self-adjoint, non-negative definite bounded linear transformation P = P* satis- 
fying (3.29). Moreover, for each e > 0, the operators [A"} l ~ l P belong to C(X,X). 

We note that Theorem 2.1 in [9] also states that c can be set equal to zero if A is self-adjoint, normal, or 
has a Riesz basis of eigenvectors. However, it is interesting to note that A is neither normal nor self-adjoint. 
For this ID problem, A does have a Riesz basis. However, this property is not needed to establish the 
following representation. 

Theorem 3.3 There exist functions Ar 2 i (’ 5 *), *22(*, *) such that 

(1) &21 ( ■> *) £ 1] * [0j 1]); ^22(*> ’) £ ^2([0 j 1] x [0 j 1]); 

( 2 ) *22(C>$) = *220,0, 

(3) the mapping t — ► &2i(00 belongs to for almost all £ € [0, 1] 
and one has the representations 


[P21<I>}{ 0 = / k 21 {C,t)d>{t)dt, <j> € Hq , 

Jo 

[P 2 2<t>m = f k 22 {C,t)<f>{t)dt, <t> £ L 2 . 

Jo 


6 



The proof of Theorem 3.3 is rather tedious and will not be included here. However, we note that the 
proof is similar to the proof of Theorem 3.2 given in [7]. It relies on a classical theorem by Fullerton (see 
Theorem 6 in [5]). However, the specific structure of the A operator combined with the special B operator, 
B = [0 Ix ] T , is needed to carry out the proof for this hyperbolic case. Theorem 3.3 leads easily to the 
following representation. 


Theorem 3.4 There exist functional gains k v ( *, ■) and &$(*,-) 

(V M-, 0 e C([0, l] X [0, i]), ■) € L 2 {[ o, i] x [o, i]) ; 

(2) k v (C,z) = k v (zX) 

(S) the mapping z — ► & 5 (£,z) belongs to H 2 for almost all ( £ [0, 1] 
and the LQR feedback control law has the representation 

u(f>C) = -j£ ks ( < ' ,z ^ w ( t ’ z ) dz - J Q MC> 

Proof: Since K : Hq x L 2 — > ► L? is given by 


«'(•) " 

(•) J 


= - f *2 

Jo 


-j: 


21 (C>z)w(z)dz- / k 22 (C,z)v(z)ds ) 


(3.35) 


let & v (-, •) = k 22 (‘, )• Integration by parts on the first integral yields (for w(-) £ Hq) 


f k 21 ((,z)w(z)dz = - [ If k 2 i(Ct)dt 
Jo Jo ]_Jo 

The representation (3.35) follows, where 


w ( {z)dz 


4,(C,z) = -/ MOO*. (3.36) 

Jo 

and z — > k s (Q } z) belongs to H 2 . 

We turn now to some numerical experiments to illustrate the representation (3.35). In addition, we 
consider other damping models and present numerical results that clearly show the role that damping plays 
in the smoothness and existence of functional gains. 


4 Numerical Experiments 

We consider the wave equation with various forms of damping. We have the representation (3.35) for Kelvin- 
Voigt damping and similar results can be obtained for structural damping. However, as we shall see below, 
viscous damping is not sufficient to ensure the existence of L 2 functional gains even though the feedback 
operator exists and is bounded. We set 7 = .25 and use standard piecewise linear finite elements to compute 
k s (-,‘) and ^(-,*). In particular, we use finite elements to approximate the Riccati equation (3.29) and 
construct K as defined by (3.30-3.32). 

In order to show convergence of the scheme we compute £^(-, -) and k^(-, •) for Kelvin- Voigt damping. 
Here, iV represents the number of elements in the model (i.e., [0,1] is partitioned into (N-fl) subintervals). 
Figure 4.1 shows the convergence of k^(-> ♦) and •) to £*(•, -) and k v (-, •), respectively. Observe that the 
convergence of Ar^(-, •) is very rapid. 

Although we presented theoretical results for the case of Kelvin- Voigt damping only, Figure 4.2 shows 
that similar results hold for structural damping. In particular, an integral representation exists and the 
finite element approximations converge. It is important to note that the functional gain for strain, •), 
remains smooth and is the same as the gain obtained with Kelvin- Voigt damping. However, there is a 
marked difference in the smoothness of k v (- } -). This functional gain has a sharper “peak” at £ = z than 
the corresponding velocity gain for Kelvin- Voigt damping. Although we expect that for structural damping, 
fc w (C> •) belongs to H l , we conjecture that t — ► k v (£, t) is not H 2 . 
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Finally, for viscous damping, the numerical results shown in Figure 4.3 show that again, (*, *) is well 
behaved and yet •) seems to be as singular measure concentrated at ( = z. Thus, we conjecture that 

the representation (3.35) does not hold for the wave equation with viscous damping for any L 2 function 
k v (- r ). Similar results were noted in [8] for undamped beam equations and in [7] for parabolic equations 
with highly unbounded input operators. 

The numerical results also indicate that the strain functional gains are independent of the damping model. 
As shown in Figure 4.4, does not change as the damping model changes. Thus, it seems from the 

numerical experiments that damping has the most impact on the existence and smoothness of k v (-, •). 

5 Conclusions 

In this paper we have provided an integral representation theorem for the LQR feedback operator for the ID 
wave equation with Kelvin- Voigt damping and fully distributed control. The motivation for this effort comes 
from the problem of designing suboptimal low order dynamic compensators and for optimal sensor location. 
Although the theoretical results presented here are restricted to problems with Kelvin-Voigt damping, the 
numerical results suggest that more general results are available. Also, these same numerical results indicate 
that once the damping is insufficient to ensure the analyticity of the control system, the existence of L 2 
functional gains is in doubt. In particular, the solution of the Riccati equation may not be Hilbert-Schmidt. 

The LQR problem for distributed parameter systems has received considerable attention during the past 
ten years. However, problems in which B , Q and R are all non-compact have not been fully explored. Temam 
[12] considers the differential Riccati equation with B = R = J, but assumes that Q is Hilbert-Schmidt. 
Similar results are found in [4, 6]. Finally, as the numerical evidence suggests, many theoretical issues are 
not yet settled. 
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k_s(.,.), N = 32 k_v(. T .), N = 32 



Figure 5.1: Convergence of the Functional Gains, Kelvin- Voigt Damping. 
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k_s(.,.), N=8, Kelvin- Voigt k_s(.,.), N^16, Kelvin- Voigt 



Figure 5.4: Structural Gains for Various Damping Types, N = 8,16. 
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